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Abstract 

Strained coherent heteroepitaxy of III-V semiconductor films such as In^Gai_^As/GaAs has potential 
for electronic and optoelectronic applications such as high density logic, quantum computing architectures, 



ON laser diodes, and other optoelectronic devices. Crystal symmetry can have a large effect on the morphol- 

ep 

ogy of these films and their spatial order. Often the formation of group IV strained heterostructures such 

as Ge deposited on Si is analyzed using analytic models based on the Asaro-Tiller-Grinfeld instability. 

However, the governing dynamics of III-V 3D heterostructure formation has different symmetry and is 

-4— > 

more anisotropic. The additional anisotropy appears in both the surface energy and the diffusivity. Here, 

+-* 

the resulting anisotropic governing dynamics are studied to linear order. The resulting possible film mor- 

£ 

phologies are compared with experimentally observed In^Gai-^As/GaAs films. Notably it is found that 

surface-energy anisotropy plays a role at least as important as surface diffusion anisotropy if not more so, 
O 



in contrast to previous suppositions. 



Strained coherent heteroepitaxy of In x Gai_ x As/GaAs at high temperatures can lead to dense 
somewhat correlated 3D film growth, for example, rolls (Figs. 2c and d in Ref. 1 !), as well as cor- 
related dense arrays of self-assembled quantum dots (Fig. 1 in Ref. 2 ) and quantum dot chains in 
multilayers (Fig. la in Ref.^. These structures are of great interest for the next generation of 
electronic and optoelectronic materials due to their quantum confinement effects, 4 and they are 
an excellent example of nanoscale self-assembly that can augment or replace traditional litho- 
graphic techniques. The density and degree of correlation in these film morphologies suggest 
that their growth is barrierless and cooperative. Furthermore, these structures form at high tem- 
perature, suggesting that the surfaces might be thermally roughened. Thus, the film morphology 
is likely governed by the Asaro-Tiller-Grinfeld (ATG) instability or related mechanism; whereby 
surface diffusion is driven by changes in elastic and surface energy. ^ 7 * 8 * 9 ! Unlike modeling group 
IV structures (e.g. Ge x S i i _ X /S i) J1QHU modeling the growth of III-V heterostructures as an ATG- 
like process requires full consideration of three contributions to anisotropic pattern formation, 
namely elastic anisotropy, surface-energy density anisotropy and diffusion anisotropy. The role 
of anisotropic effects has been increasingly recognized in nanoscale epitaxial self-assemly ] 12 ! 13 ! 1 ^ 
Here, a linear model is presented that has appropriate symmetries for III-V structures generally, 
and In x Gai_ x As/GaAs in particular. 

Similar models have been largely applied to group IV heteroepitaxy, such as Ge x Sii_ x /Si, that 
have four-fold rotational symmetry of the crystal surface, but III-V systems have only two-fold 
rotational symmetry. The experimental observations bear out this difference. In x Gai_ x As/GaAs 
structures grown on (001) surfaces form rows or rolls that are aligned closer to the [1 10] direction 
than the [110] direction. On the other hand, Ge x Sii_ x /Si structures form a four-fold symmetric 
pattern with alignment along the (100) directions. 15 This spatial alignment of In x Gai_ x As/GaAs 
structures closer to the [110] direction has been attributed to differences in surface diffusivity^EE] 
and possibly to packing effects and surface faceting.^ To test these hypotheses, one must obtain 
similar structures from a complete model of surface evolution. A linear analysis of surface dif- 
fusion based on the original ATG instability has been used to explain anisotropic patterning in 
Ge x Sii_ x /Si heteroepitaxy] 1 °l 1 1 1, and it would be remiss not to apply it to In x Gai_ x As/GaAs. In 
group IV systems, the symmetry requirements forbid anisotropic surface-energy and diffusion ef- 
fects to linear orderpSI Thus, for Ge x Sii_ x /Si the (100) alignment comes from elastic anisotropy. 
The lower symmetry of III-V (001) surfaces allows surface-energy and diffusion anisotropy ef- 
fects; thus, all three contributions interact to determine the alignments of 3D morphological fea- 



tures. 

In the following, a continuum linear model of anisotropic film evolution is presented including 
the general framework, then energetic and diffusional sources of anisotropy. Then the conse- 
quences for film morphology are presented and compared with experimental observations. 

Early stage small fluctuations in film height are an important determiner for final film morphol- 
ogy ^ and can be analyzed using linear dynamics of the film height that is a sum of the average 
film height (JF) and film height fluctuations (h(x,t)), Jf?(x,t) = J^ + h(x,t). The normal modes 
of the film evolution are periodic; thus, Fourier components of the height fluctuations are used; 
h k =A~ l f d 2 xe~' k ' x h(x) and h (x) = £k e ' k x h k , where A is the area of the substrate, and periodic 
boundary conditions are assumed. Following Ref. 17 , each Fourier component is initially governed 
by energy dissipation and thermal fluctuations, 

d t h k = o(k)h k + ik-y/2k b TT>-ri k , 
CT(k) = -(k-D-k)/(k), (1) 

where <r(k) is the dispersion relation, D is the surface diffusivity, r} k (t) is a normally distributed 
stochastic variable with mean (%(?)) = 0, and variance (^k(?)^(?')) = A~ l d^wS (t — t'). 5kk' 
is the Kronecker delta, 8(t — t') is the Dirac delta, and /(k) is the energy cost function, formally 
defined as the second derivative of the free energy per unit area {& /A) with respect to each Fourier 
component, 

f(k) = d hk d ht (^/A) (2) 

From <r(k), one can determine various length and time scales as well as pattern orientation and 
alignment. 

The film free energy cost function consists of two terms due to elastic strain and surface energy 
/(k) = /eiast.(k) +/surf. (k). For now, the wetting energy contibution 18 is neglected for simplicity. 
For an (001) surface, elastic anisotropy leads to 4-fold-symmetric surface energetics and dynamics, 
and surface energy anisotropy leads to 2-fold symmetric energetics and dynamics as explained 
below. 

The effect of elastic anisotropy alone has been discussed previously] 10 * 1 1 I 16 I S 1 an d the most 
important results are summarized below. For (001) surfaces, the elastic part of the energy cost 
function can be approximated as-^ 

/elast(k) = - [<?o° cos 2 (2%) + ^5° sin 2 (20k)] k, (3) 



where <§§° and ^45° are constants, 0k is the angle between k and the [100] direction, and k = ||k||. 
/elast. is most negative along the (100) directions; thus ripples perpendicular to these directions 
release the most elastic energy. Without competing anisotropic effects, the anisotropy of / e i ast , 
causes initial alignment of 3D gratings along the (100) directions] LEUI! p or i n As, $0° = 8.13 x 
10 9 erg/cm 3 , and <%° = 6.92 x 10 9 erg/cm 3 . 

One can derive the general form of / sur f as follows. Assume that the total free energy is an 
integral in the x — plane over a local free energy kernel ft) that depends on the surface orientation 
VJf , J^urf. = fd 2 xco(VJ4?(x)) . In terms of the surface energy density 7, ft) = y(Vjf )[1 + 
(VJf (x)) ] 1 / 2 . Using Eq. UA the free energy cost function is 



/ surf .(k) = k-(»"-k, (4) 

where a>" = <?v^tt)|v^=o — 7(0)1+7" (0), co" acts like a surface stiffness stabilizing short wave- 
lengthsP 

From Eq. |4j one can determine the possible form of any anisotropy. a>" is a rank 2, dimension 
2, symmetric tensor, and it has two eigenvalues and eigenvectors. Due to crystal symmetry, namely 
reflections through the (110) and (l 10) planes, the eigenvectors of co" must be in the [110] and 
[110] directions with eigenvalues, ft)nio] an ^ Wfiioi • The surface energy anisotropy coefficient 
can be defined as 5 surf . = ft> [ll0] /ft> [110 ] so that a>" = ooq (n [no] n [no] + 5 surf n [ll0] n [ll0] J, where 
COq = 0)[no], nil is the unit vector in the specified direction, and the vector products are outer 
products. Thus, 5 sur f = 1 indicates an isotropic surface stiffness, and 5 sur f. > 1 (< 1) indicates 
an enhanced stiffness in the ±[110] (±[110]) directions. In terms of the surface energy density, 
5 S urf. = (7(0) + 7[iIo])/(7(0) + Tfno]), where f[...] are the eigenvalues of f. 

The exact values of COq and 5 sur f. are not known, so OJo will be taken to be about 770 erg/cm 2 .123 
Then, taking in turn 5 sm -f = 1, 1.5, 2, and 5 the total energy cost functions are shown on the left- 
hand column of Fig. [TJ Only values of 5 sur f. > 1 are shown, corresponding to enhanced stiffness 
in the ±[110] directions. 

In all cases, the energy cost function has negative minima indicating the wavevectors of the 
energetically most unstable modes, k#. In agreement with previous findings LMUJ for the isotropic 
surface energy case, k^ point along the (100) directions with magnitude Ice = <£b°/ (2ft>o) . How- 
ever, for <5 sur f. > 1, the peak wave vectors move towards the ±[1 10] directions, and the magnitudes 
decrease so that kg < <%V (2&>o). For <5 sur f. > 2.4, the peaks actually merge to form one peak with 
k^ along the ±[1 10]. An opposite trend would be found for 5 sur f. < 1 (not shown). 
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Figure 1 : Graphical "multiplication table " showing effects of diffusional and surface energy anisotropy. 
Top row shows (k-D-k)/Do for various surface-diffusion anisotropics cW.. Left column shows /(k) for 
various surface-stiffness anisotropics 5 sur f.. Interior grid shows resulting dispersion relation a(k) (Eq. [I]). 
Horizontal direction (— ►) is [100]. Vertical direction (f) is [010]. All axes span ±lrad/nm. &s indicate k^ 
for /(k) and ko for o~(k) plots. Dashed lines perpedicular to ko indicate mean alignment of grating rows or 
rolls (Fig. |2). 



In addition to energetics, one must consider dynamics. The surface diffusivity is a rank 2, 
dimension 2 tensor, and it must have the same symmetry as the crystal surface. Thus, it has 
two eigenvectors and eigenvalues corresponding to the [1 10] and [1 10] directions just as &>" does. 
Let S di f. = D [110] /D [ll0] and D = D [uo] . Then, D = Do(5^n [110 ]ii [110] +n [ll0] n [ll0] ). 5 dif . > 
1 indicates slower diffusion in the [110] direction, which is believed for As terminated (001)5 
InAsPThe top row of Fig. Fjl shows (k-D-k) /D for 5 dif . = 1, 2, 5, 10 and 100. 
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Figure 2: "Simulated" film heights for various values of 5dif. and 5 sur f. corresponding to a(k)'s appearing 
in Fig. [I] All figures oriented with [100] horizontal (— =>) and [010] vertical (|). Dashed lines show mean 
grating or roll alignment predicted by a(k). 

Three contributions to anisotropy in surface morphology have been discussed. The first two, 
elastic anisotropy and surface-stiffness anisotropy contribute to an anisotropic energy cost func- 
tion, /(k). The last, diffusional anisotropy, multiplies /(k) to give an anisotropic dispersion rela- 
tion <j(k) (Eq.[T]). The peaks in <j(k) give the fastest growing perturbation wavevectors, krj. Fig.[T] 
shows a figurative multiplication table of the combined effects of surface-stiffness anisotropy and 
diffusive anisotropy on the dispersion relation. The left column and top row show the surface- 
stiffness and diffusivity anisotropy respectively while the inner grid shows the combined effects. 
Resulting real-space morphological alignements for each case will perpendicular to ko and are 
shown by dashed lines in Fig. [I] 

From examination of Fig.[T]and calculations using other <5 sur f. and 8^. values (not shown), one 
finds: 



1. When only elasticity is anisotropic (<5 sur f. = <5dif. = 1), the peaks are oriented along the 



(100) directions, the same as the peaks in /(k), but the magnitudes &o = |ko I = (3/2) k£ = 



3c0q/ (4<f °) as expectedP 21 

2. When 5<jif, > 1 and 5 sur f. = 1, diffusivity is greater in the "fast" [110] direction, but with 
isotropic surface stiffness the peaks move towards the fast ±[110] directions, and the mag- 
nitude ko decrease. 

3. When <5dif. = 1, but <5 sur f. > 1, the peaks move towards the "slow" ±[110] directions. In 
particular, for <5 sur f. > 1.94, the pairs of peaks in the ±[110] directions each merge into one 
peak exactly along the ±[110] directions. Note that the energetic minima at k £ merge at 
values, <5 sur f. > 2.4. 

4. When 5dif. > 1, and 5 sur f. > 1, the two effects compete. For example, for 5dif. = 100 and 
<5 SU rf. = 2, the peak positions ko appear almost exactly along the the (100) axes, and for large 
<5surf. an d <5dif., for example, <5 sur f. = 5 and <5dif. = 100, the merged peaks in /(k) are split by 
the anisotropic surface diffusivity back into four peaks that appear near the slow ±[110] 
directions. 

The initial film morphology will be quasiperiodic with reciprocal lattice vectors given by the 
peaks in c(k) (Fig. [I]). For visualization, one can sample film height fluctuations from the ap- 
propriate statistical distribution. Following RefP-^, Eq. [T]is used to find the ensemble means and 
variances of the individual film height Fourier components. Taking the ensemble average of Eq.[T} 
dt (hk) = c(k) (/%)• Since initially (h^) = 0, (/?k) = for all t. Using the stochastic chain rule 
for d t \hk\ and taking an ensemble average, d t ( |%| ) = 2<7k ( \hk\ / + 1k\,T (k • D ■ k) /A. Solv- 



ing, (|/?k| 2 ) = hT ( e 2<T ( k )W — l j /[A/(k)], where ^ge is chosen so that the r.m.s film height 
fluctuation is an appropriate value for small fluctuations, /? r . m .s = (Lk(l^k| 2 )) = 2.5 A. 

Fig. [2] shows sampled initial film morphologies corresponding to the dispersion relations c(k) 
from Fig. [T] Dashed lines in both figures are perpendicular to ko and show the mean orientations 
of either grating rows or rolls. Dispersion relations (c(k)) with four peaks lead to quasiperiodic 
gratings that are suitable precursors to self-assembled quantum dots. Dispersion relations with 
two peaks lead to roll structures. These structures, however, might evolve significantly during 
later development perhaps in a fashion similar to ripening. Such evolution might depend in detail 
on the surface energy density y(V^f ) to higher than linear order and is thus a subject for future 
investigation. 

7 



Figs, [l] and |2| indicate that fast diffusion in the [110] direction alone can not be responsible 
for the observation of film morphologies aligned either along or close to the [110] direction] LJJ2]3] 
In fact, it gives the opposite result, alignment close to the slow [110] direction. [110]-aligned 
patterns occur either for moderate values of 5 sur f. and small or zero 5dif., or for the case where 
both 5 sur f. and 5dif. are large. It is generally accepted that the [110] direction is the fast diffusion 
direction 21 suggesting that the latter case is more likely than the former. As mentioned above, 
<5 S urf. < 1 would have the opposite effect; thus observations of SAQD arrangements can be used to 
infer that 5 sur f > 1. A Fourier transform of Fig. 1 from Ref. 2 has broad peaks that overlap the slow 
[110] directions significantly. Perhaps the spectral support along the [110] direction is initially 
suppressed by diffusional anisotropy and appears later due to non-linear effects which will be the 
subject of future investigations. The other mentioned images from Refs. 1 3 have spectral support 
that is close to the [1 10] direction but with less overlap. 

A linear surface evolution model was presented that incorporates elastic, surface energy and dif- 
fusional anisotropy. It was found that fast diffusion in the [1 10] direction gives initial morhpology 
aligned along the slow [110] direction. Thus, oberved [110] directed alinement of morphology^EE] 
must be caused in whole or in part by surface energy anisotropy. Future non-linear modeling is 
needed to confirm this observation and may reveal even more complicated mechanisms. 
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